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ABSTRACT 

We examine the left -right asymmetry in the cleaned COBE/DIRBE near- infrared 
data of the inner Galaxy and show (i) that the Galactic bar is probably not seen very 
nearly end-on, and (ii) that even if it is, it is not highly elongated. The assumption of 
constant mass-to-light ratio is used to derive simulated terminal-velocity plots for the 
ISM from our model luminosity distributions. By comparing these plots with observed 
terminal velocities we determine the mass-to-light ratio of the near-IR bulge and disk. 

Assuming that all this mass contributes to gravitational microlensing we compute 
optical depths r for microlensing in Galactic-centre fields. For three models with bar 
major axis between 10° — 25° from the Sun-Galactic Center line, the resulting optical 
depths in Baade's window lie in the range 0.83 x 10~^ ^ t < 0.89 x 10~^ for main- 
sequence stars and 1.2 x 10~^ ^ t < 1.3 x 10~^ for red-clump giants. We discuss 
a number of uncertainties including possible variations of the near-infrared mass-to- 
light ratio. We conclude that, although the values predicted from analyzing the COBE 
and gas velocity data are inconsistent at the 2 — 2.5a level with recent observational 
determinations of r, we believe they should be taken seriously. 
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1 INTRODUCTION 

It is widely recognized that observed optical depths to mi- 
crolensing along various lines of sight from the Sun con- 
stitute important constraints on models of the Milky Way. 
In his seminal study of microlensing Paczynski (1986) es- 
timated optical depths r for lensing under the assumption 
that the Galaxy is axisymmetric. He found that in fields to- 
wards the Galactic centre r is ~ 10~®, and numerous stud- 
ies have since confirmed the correctness of this estimate. 
When T was actually measured, significantly larger values 
have been reported - r = (3.3 ± 1.2) x 10"® (la) by (Udal- 
ski et al (1994) and r = (1.9 ± 0.4) x 10"® by Alcock et 
al. (1995)]. It was suggested (Paczynski et al. 1994; Zhao, 
Spergel & Rich 1995) that the large measured values of r 
reflect the fact that the Galaxy is barred (see, e.g., reviews 
in Blitz & Teuben 1996). 

In principle a suitably oriented and elongated bar can 
significantly enhance r. To see this, imagine deforming an 
initially axisymmetric bulge into a bar that points nearly 
to the Sun in such a way as to hold constant the velocity 
dispersion within Baade's window [{I, b) = (1°, —3.9°)]. This 
deformation will increase both the surface density of stars 
towards Baade's window and the depth along this line of 
sight within which the stellar density is high. Since a star 
offers the largest cross section for lensing a background star 
when it lies halfway between the background star and us. 



both of these factors will increase r. However, in this paper 
we show that if one assumes that near-infrared luminosity 
density is a fair tracer of stellar mass density in the inner 
galaxy, then r is rather precisely constrained to be r = 10"® 
by a combination of (i) the near and far-infrared surface 
brightnesses that were measured by the DIRBE experiment 
aboard the COBE satellite, and (ii) the kinematics of gas in 
the inner few kiloparsecs as reflected in longitude-velocity 
plots from radio-frequency emission-line studies of the ISM. 

In substantiating our claim we rely heavily on results 
obtained in three previous papers, namely Binney & Ger- 
hard (1996; Paper I), Spergel, Malhotra & Blitz (1997; Pa- 
per II) and Binney, Gerhard & Spergel (1997; Paper HI). 
Paper II cleans the COBE/DIRBE data for the effects of 
obscuration by dust, while Papers I and HI develop and ap- 
ply a Lucy-Richardson algorithm for recovering models of 
the Galaxy's three-dimensional luminosity density from the 
cleaned near-IR surface photometry. 

This paper is organized as follows. In Section 2 we fo- 
cus on left-right asymmetry within the observed near-IR 
surface brightness measurements as a key diagnostic of the 
morphology of the Galactic bar. In particular we explain 
how the length, axis ratios and orientation of a bar are re- 
flected in its distribution of left-right asymmetry. We then 
compare these asymmetries both with the asymmetry that 
is apparent in the COBE data, and the asymmetries that 
are predicted by our models of the Galaxy. These compar- 
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isons suggest that the Galactic bar is not seen very nearly 
end-on, and even if it is, it is not highly elongated. 

In Section 3 we assume that our models of the near- 
IR luminosity density can be converted into mass models 
by multiplying by an appropriate mass-to-light ratio T. We 
calibrate T by comparing simulated longitude-velocity plots 
of the ISM with observed {l,v) plots. Finally we determine 
T. In Section 4 we discuss uncertainties and the possible 
causes for the difference between the optical depth inferred 
here and that inferred from the microlensing experiments. 



2 LEFT-RIGHT ASYMMETRY 

As Blitz & Spcrgcl (1991) have emphasized, the key to using 
photometry to detect departures from axisymmetry within 
the Galaxy is the study of the ratio A(Z, 6) that is defined 
for positive I by 



A{l,b) = I{l,b)/I{-l,b), 



(1) 



where I{l,b) is the surface brightness of the Galaxy at lon- 
gitude I and latitude b. Fig. 1 is a contour plot of A for 
the L-band COBE data after the latter has been corrected 
for dust obscuration as described in Spergel, Malhotra & 
Blitz (1996). Values of A are available for pixels that are 
1.5° square, but the values contoured in Fig. 1 are those 
obtained by smoothing the raw data by 5 pixels in I and 2 
pixels in b with the SMOOFT routine of Press et al. (1986). 
We shall call a plot such as Fig. 1 an asymmetry map. 

To appreciate the significance of Fig. 1 it is instructive 
to construct analogous figures for simple Galaxy models. 
Fig. 2 shows several illustrative cases for the model intro- 
duced in Paper I. This consists of a bar/bulge and a disk. 
Its luminosity density is given by 

(2a) 



jix) =io[/b(x)H-/d(x)], 



where 



Bo 



al^vC (l + o/ao)i « 
/d = ( ha 



zo 



2 V Z 

^ + ^ + 72 
^2 1^2 



2\ 1/2 



and R={x^+yy/''. 



(2b) 



There are four important parameters associated with this 

model bar: the length am, the axis ratios jy, and the angle 
^0 between the Sun-centre line and the long axis of the bar 
(which is defined such that the near end of the bar lies at 
Z > for < 00 < 90°). The other parameters may be 
held constant at the values given in Table 1. We assume 
throughout that the Sun lies 14 pc below the Galactic plane 
as was deduced in Paper III. 

The top-left panel of Fig. 2 shows the asymmetry map 
for fflm = 1-91 kpc, r] = 0.5, C, = 0.3, 4>o = 25°. In contrast to 
the data map. Fig. 1, all contours are positive, that is, for 
the given parameters, the model (2) is nowhere more than 
0.05 mag brighter at Z > than it is at Z < 0. The brightness 
contrast peaks near {l,b) = (8°, ±4°). For the middle-left 
panel of Fig. 2 the vertical axis ratio C has been doubled to 
C = 0.6. This moves the peaks in the asymmetry map out 
away from 6 = and increases their intensity. The top-right 



Figure 1. The asymmetry map of the COBE L-band data. Con- 
tours arc spaced by 0.05 magnitudes. Dotted contours indicate 
that the Galaxy is brigiiter on the right than on the left. 

Table 1. Parameters of the models 
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panel in Fig. 2 shows the effect of decreasing the planar axis 
ratio from 0.5 to 7? = 0.3 whilst holding the vertical axis 
ratio constant at (" = 0.3. This leaves the peaks at approx- 
imately the same location but makes them more intense. 
The middle-right panel of Fig. 2 shows the effect of dou- 
bling the length of the bar from 1.91 kpc to 3.82 kpc whilst 
holding the axis ratios constant at 0.3. This both moves 
the peaks towards larger I and higher b and makes the sur- 
rounding contours significantly more extensive. Comparison 
of the middle and bottom-left panels of Fig. 2 shows the ef- 
fect rotating the bar to a more nearly end-on orientation; the 
bottom-left panel is for </>o = 10° rather than 00 = 25° . This 
leaves the peaks in approximately the same location because 
even at 00 = 25° the line of sight through the Galactic cen- 
tre effectively passes down the length of this fairly fat bar 
rather than intersecting it transversely. In this regime the 
magnitude of the peaks decreases with 4>o - by symmetry it 
must vanish at 0o = 0. 

The bottom-right panel of Fig. 2 shows the asymmetry 
map of a long (am = 3.82 kpc) thin (77 = 0.25, C, = 0.2) bar 
seen nearly end-on (0o = 10°). Comparison of this asymme- 
try map with that shown in the middle-right panel of Fig. 2 
reveals that the effect of rotating a long thin bar to a more 
nearly end-on orientation is two-fold: it decreases the mag- 
nitude of the peaks as in the case of a short bar, but it now 
also shifts them to slightly lower longitudes. This is because 
when a thin bar is oriented at 0o = 25°, the line of sight 
through the Galactic centre intersects the bar transversely. 
Finally, comparison of the two lower-right panels with the 
corresponding lower-left panels shows that rotation of a bar 
towards end-on orientation combined with elongation of the 
bar can leave the peaks in the asymmetry map at approx- 
imately the same locations. However this combination of 
rotation and elongation modifies the shapes of the contours 
that surround the peaks in a characteristic way: they be- 
come more elongated in the longitude direction, and they 
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Figure 2. Asymmetry maps of analytical models. The parameters in eqs (2) of each model are given at the top of its panel. The 
contours are as in Fig. 1. 



Figure 3. Asymmetry maps of iterated models. The parameters in eqs (2) of the model from which the iterations started are given at 
the top of each panel. The contours are as in Fig. 1. 



Figure 4. Projections parallel to z of three final models: From left to right: for = 25° started from a short-fat analytic bar (middle 
left panel of Fig. 2), and started from a longer bar of similar cross-section (middle-right panel of Fig. 2); for <^o = 10 started from a long 
thin bar (lower right panel of Fig. 2) . 



turn upward towards high latitudes at I > 20°. This effect 
arises because the far end of a long near end-on bar appears 
small because it is far off, so away from the plane and at 
large longitudes there is no counterpoint to the brightness 
contributed by the apparently large near end of the bax. 
Conversely, at small I and h light from the far end of the bar 
provides an effective counterpoint to light from the near end 
of the bar, and the asymmetry is small. 

Fig. 3 shows the asymmetry maps of the models that 
are generated by the algorithm of Papers I and 111 from 
the COBE Z/-band data when the Lucy-Richardson itera- 
tions are started from the analytic models that generate the 
asymmetry maps of the corresponding panels in Fig. 2. Sev- 
eral points are noteworthy in this figure: 

(i) Especially in the 0o = 25° case, the asymmetry map 
of the model generated by the algorithm is remarkably inde- 
pendent of the initial analytic model, even though the length 

and axis ratios of the latter vary by factors of two. 

(ii) The 0o = 25° final model provides a remarkably 
good fit to the global morphology of the COBE asym- 
metry map of Fig. 1. In particular, in the middle-right 
panel of Fig. 3 the two long ridges of mgiximal asymmetry 
around (Z, 6) = (10°, ±6°) strongly resemble ridges in Fig. 1 
that have similar locations and orientations. Moreover, the 
height of these ridges is nearly the same in the two figures, as 
is the general shape of the low-lying contours near I = 20°. 
Finer details of the observed asymmetry map shown in Fig. 1 
are not reproduced in Fig. 3. These include the sharp local 
maximum at {l,b) — (16°,— 1°) and the strong north-south 
asymmetry. The former is almost certainly associated with a 
feature in the disk and does not concern us here. The north- 
south asymmetry in the observed asymmetry map could be 
a sign that the bar is slightly inclined with respect to the 
plane, as radio-frequency emission lines from the ISM sug- 
gest (e.g., Liszt & Burton 1996). 

(iii) The asymmetry maps that are shown in the bottom 
panels of Fig. 3 for the case (t>u = 10° both deviate from 
the observed asymmetry map of Fig. 1 in having contours 
at {l,b) ~ (25°. ±7°) which slope towards higher I and |6| 
rather than towards lower / at higher |6|. Thus whereas in 
the case of 4>o = 25° the Lucy-Richardson iterations are able 
to eliminate the extensions towards high I and |&| that are 
apparent in the middle-right panel of Fig. 2, the iterations 
cannot eliminate this unsatisfactory feature of the initial 
models in the case 4>o = 10°. When a bar is viewed nearly 
end-on it is inevitable that the sky is brighter up and to the 
left than it is at the corresponding point on the right. 



3 THREE-DIMENSIONAL STRUCTURES 
AND MASS DENSITIES 

Fig. 4 shows projections parallel to z of the bars that were 
generated by six Lucy-Richardson iterations on the COBE 
L-band data from the analytic bars whose asymmetry maps 
are shown in the middle- left, middle-right and bottom- right 
panels of Fig. 2. These bars are very similar to one an- 
other. In particular, the 25°-model that started with a long 
bar finishes (middle panel) with a bar of exactly the same 
length as the 25°-model that started with a bar of half the 
length (left-hand panel). Comparing the middle two panels 
of Figs 2 and 3, it is clear that the Lucy-Richardson iter- 
ations have dealt with the problem posed by the extension 
towards large I and |6| of the peaks in the middle- right panel 
of Fig. 2 by shortening the bar. The 10°-model shown in the 
right-hand panel of Fig. 4 has a slightly shorter bar than do 
the 25°-models shown in the middle and left-hand panels. 
Thus again the Lucy-Richardson iterations have shortened 
the bar, but the bottom-right panel of Fig. 3 shows that 
now shortening is not accompanied by elimination of the 
unwanted extensions of the asymmetry peaks towards high 
I and |6|. On the other hand, shortening the bar has made it 
possible to broaden it, and altogether this has increased the 
amplitudes of the peaks (cf. Fig. 2), as the observations de- 
mand. These experiments demonstrate that the gross struc- 
ture of the models recovered by Lucy-Richardson iteration 
can be understood in simple general terms. Hence we may 
confidently assert that the Galactic surface-brightness dis- 
tribution is such that the length of the bar that is required 
to fit the COBE data for given (f>o is nearly independent of 
<po for (f> < 30° . We show below that the same is true of the 
optical depth to microlensing. 

Before an optical depth can be determined, it is neces- 
sary to associate a mass-to-light ratio T with a model. This 
we do as follows. We first construct a luminosity model for 
the entire Galaxy out to the solar radius, by combining the 
luminosity distribution recovered by the Richardson Lucy 
algorithm inside our 5kpc x 5kpc x 1.4 kpc computational 
box with the initial model of Paper III in the region outside 
the box. Then we reconstruct the cusp in the centre, which 
is unresolved in the DIRBE data, by fitting a power law to 
the multipole expansion of the luminosity distribution in the 
range 350 — 500 pc. We next evaluate the combined model's 
potential under the assumption of constant T. Then we use 
an SPH code to simulate the fiow of gas inside ~ 8 kpc [see 
Englmaier & Gerhard (1997) for details], including a higher 
resolution simulation for the central ~ 4 kpc. In this simular 
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tion, the assumed pattern speed of the bar is 60 kms-Vkpc, 
placing corotation at 3.1 kpc. The gas flow velocities in the 
simulation are transformed into a putative local standard of 
rest (LSR) frame at the position of the Sun at i?o = 8 kpc. 
For this transformation we assume that the LSR has a tan- 
gential velocity vq, but no component of motion in the direc- 
tion of the Galactic center. For each line-of-sight we evaluate 
the maximum velocity that would be seen in radio- frequency 
emission-line measurements. 

We compare the resulting model terminal- velocity curve 
with the observed HI and CO terminal velocities (corrected 
for line width, but not for peculiar LSR motion; Bur- 
ton & Liszt 1993, Clemens 1985). Specifically, for vq = 
180, 190, 200, 210, 220 km s"^ we choose the value of T that 
gives the best eye-ball fit to the observed velocities in the 
region between I = 17 — 48°. Finally we choose an optimal 
(vo,T) combination. 

The thin full curve in Fig. 5 shows the fit between model 
and data for a model with (f)o = 20° under the assumption 
vo = 200 km s"^ - the fit to the data at (17° < I < 48°) 
is consistent with the scatter in the data points (~ 5%). 
In the range (9° < / < 17°) the model curve falls below 
the data by < 20kms~^. This may reflect that fact that 
the photometric inversions show a strongly elliptic disk at 
(1.5 kpc < r < 3.5 kpc), and, as discussed in Paper III, 
this could be caused by the presence of supergiants in star- 
forming regions along prominent spiral arms. If so, our as- 
sumption of constant T may not be correct in this region 
and the gravitational forces in our models may be insuffi- 
ciently accurate at radii corresponding to (9° < / < 17°). 
At / ~ 3° the thin, full model curve in Fig. 5 falls below the 
sharply peaked data. This may in part be due to the model's 
pattern speed being not quite correct, but in large measure 
it undoubtedly reflects the limited spatial resolution of our 
models: the thick, full curve in Fig. 5 shows the predictions 
for exactly the same model as the thin, full curve, but cal- 
culated from a higher-resolution SPH simulation. It can be 
seen that the thick curve halves the shortfall of the thin 
curve with respect to the data. However, this gas flow is 
not stationary, and as more gas falls towards the Galactic 
center and the vicinity of the cusped orbit is depopulated, 
the maximum velocities seen in the bulge region decrease 
again. 

The dotted curve in Fig. 5 shows the predictions of the 
model with 0o = 10°. Again vo = 200kms~^ has been 
assumed. The match to the data for this model is reason- 
able but somewhat inferior to that obtained with the corre- 
sponding 00 = 20° model (thin, full curve). Specifically, for 
(f)Q = 10° there are systematic deviations ~ lOkms"^ from 
the observed velocities. About half the difference between 
the terminal velocity curves of the 10° and 20° models is due 
to the difference in viewing angles; the remaining discrep- 
ancy is due to intrinsic differences in the model gas fiows. 

The quality of the fit between the thick, full curve and 
the data in Fig. 5 indicates that to first order we may assume 
that the mass-to-light ratios of the bulge and disk are equal. 
With this assumption we find 

p{x) = TLjL{x), (3a) 
where jL{x) is in COBE/DIRBE luminosity units per kpc^ 



Figure 5. Comparison of observed and model terminal veloci- 
ties. The diamonds show the terminal velocities and their error 
bars from the CO data of Clemens (1985). The squares show 
unpublished HI terminal velocity measurements of Burton. The 
stars are HI terminal velocities read from the figures in Burton 
& Liszt (1993). The curves show the predictions from three SPH 
gas fiow models after subtracting the component of solar tangen- 
tial motion, DgsinZ, with vq = 200 kms"'^. Thin full curve: gas 
flow in the mass distribution obtained by inverting the COBE 
surface brightness map with 0o = 20°, at intermediate numerical 
resolution. Thick full curve: the same model in a higher resolu- 
tion simulation. Dotted curve: gas fiow in model inverted with 
<f>0 = 10°. The Sun's galactocentric radius has been taJcen to be 
Ro = 8 kpc. The models have been observed at time t = 0.12 Gyr 
and have been scaled to fit by eye the observed data between 
I = 17-48°. 

as in Paper III and 

Ti(0o = 20°) = 3.47 X 10* M©/COBE unit, 

(3b) 

rL{(f>o = 10°) = 3.37 X 10® M0/COBE unit. 

Thus the derived mass-to-light ratio Ti is almost indepen- 
dent of (po- 

The uncertainties in the normalization of p{x) arc as 
follows: Changing the tangential velocity of the LSR to 
180 kms"^ and 220 kms"^ changes Tl by -10% and -|-10%, 
respectively. From the above discussion of the velocity peak 
in Fig. 5, we conclude that the near-IR mass-to-light ratio of 
the bulge can differ from that of the disk by at most < 20%. 
Other uncertainties such as that caused by the unknown 
projection angle 0o are smaller. 

In Table 2 we give approximate bulge masses interior to 
cylindrical radius R = 2.4 kpc for the three COBE models 
of Fig. 4. These were estimated from the complete density 
field p{x) in the following way. We associated with the disk 
all mass that lies within 0.1 radian of the plane (as seen 
from the Galactic centre). For each bounding radius R we 
compared the total mass MtotiR) inside R with the mass 
Mb(il) inside R that is, by the above definition, non-disk 
mass. We found that Mb saturates for R between 2 — 3 kpc. 
For R < 2.5 kpc the difference Mtot — Mb is approximately 
equal to the mass of the initial analytic disk model of Paper 
HI that lies within R, while further out there is significant 
disk mass that is not accounted for by the initial analytic 
model. From these results we infer that a reasonable esti- 
mate of the bulge's mass is given by the difference between 
Mtot(2.4kpc) = 1.9 X 10^° Mq and the mass of the initial 
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disk that lies at i? < 2.4 kpc. This is the quantity that is 
given for each model in Table 2. 



4 MICROLENSING OPTICAL DEPTHS 

The microlcnsing optical depth t(Ds) for a source at 
distance Ds is the probability for it to fall within one Ein- 
stein radius of any intervening star: 



10 

4nG . 



' 4nGp{Dd) DiDj 
c2 A 



dDd 



Da / p(x)x(l — x)6.x, 



(4) 



where p{Di) is the mass density of lenses at distance D^, 
X — Dd/ Ds, and Ds = Dd+Dds- t depends only on the mass 
density of lenses, and is independent of the lensing mass 
spectrum. The measured optical depth in the direction [l, b) 
is the function t{Ds) averaged over all observable sources 
that are brighter than some apparent magnitude m in a cone 
of small 51, Sb around {l,b). Following Kiraga & Paczynski 
(1994) it is customary to parameterize the distribution in 
Ds of these sources as 



= constant x p{Ds)d!+^''. 



(5) 



This assumes that the luminosity function of sources is con- 
stant along the lino-of-sight and that the fraction of stars 
brighter than luminosity L is proportional to L'^; the factor 
Ds accounts for the increase of volume with distance. Thus 



jPr{Ds)p{Ds)D!+'^ dDs 
j;°p{Ds)D's+"' dDs 



(6) 



Table 2 lists optical depths r_(j = (r)/10~^ in Baado's 
window [{l,b) — (1°,— 3.9°)], for the three luminosity mod- 
els that are shown in Figs. 2-4. In all cases the integrals 
were evaluated for the density p(i?s) = '^LjhiDs) with the 
scaling of equations (3). Two values of (3 were considered: 
P = Q, appropriate for stars that can be detected at any dis- 
tance along the line-of-sight, such as the clump giant stars, 
and (3 = —1, which is more appropriate for main-sequence 
stars (Kiraga & Paczynski 1994, Zhao, Rich & Spergel 1996, 
Zhao & Mao 1996). Figure 6 shows optical depths along the 
minor axis I = 0. The differences between models arc every- 
where small and can be quantified in first order by Eissigning 
a multiplicative constant to each model. The difference be- 
tween the /? = and P = —1 cases is of order 35% in Baado's 
window. Figure 7 shows an optical depth map for one of the 
models. Maps for other models have similar structure. 

In comparing these results with observations note that 
the self-lensing contribution of the Galactic disk is auto- 
matically taken into account. If we restrict the source stars 
to be bulge stars, the optical depth increases by of order 
25%. For a bulge star at 8 kpc distance in Baade's window, 
foreground lenses with galactocentric distances — 2 kpc, 
2 — 3 kpc, 3 — 4 kpc, 4 — 5 kpc and 5 — 8 kpc contribute 
T_6 = 0.42,0.19,0.15,0.12,0.21 to the total optical depth 

T-6 = 1.1. 

The main results from Table 2 and Figs. 6-7 are the 
following: 

(i) The predicted values of r are almost independent of 
the orientation 0o of the bar and of the initial model from 



Table 2. Microlensing optical depths in Baade's 
window for the three COBE models shown in 
Fig. 6, in units of IQ-^. 



Model 


00 


am 


M6,2.4 


«=() 


d=-l 


1 


25° 


1.91 kpc 


8.6 X 10'' 


1.2 


0.86 


2 


25° 


.3.82 kpc 


8.1 X 10" 


1.2 


0.83 


3 


lo- 


:i.82kpc 


7.2 X 10'' 


1.3 


0.86 



Figure 6. Microlensing optical depth for several models along 
the minor axis {I = 0). The short-dashed lines show optical 
depth as a function of latitude for the (t>o = 25° -bulge-disk model 

(top-left panel of Fig. 3). The long-dashed lines show optical 
depths for the long 4>q = 25° -model (middle-right panel of Fig. 3 
and middle panel of Fig. 4). The solid lines show optical depths 
for the 00 = 10°-model (bottom-right panel of Fig. 3 and right 
panel of Fig. 4). For all models, the upper curve is for (3 = and 
the lower curve is for f} = —1. 



Figure 7. Microlensing optical depth map for the model shown 
in the middle-right panel of Fig. 3 and middle panel of Fig. 4. 
The thick contour has T_6 = 1; the spacing between contours 
corresponds to a factor of 1.5. 

which the iterations were started. 

(ii) For /3 = — 1, the predicted optical depth in Baade's 
window is r_6 = 0.8 — 0.9; for comparison, the inferred 
OGLE and MACHO values for main-sequence stars are 
r-6 = 3.3 ±2.4 (2cr) from 9 events (Udalski et al. 1994), and 
r-6 = 1.9±0.8 (2cr) from 41 events (Alcock et al. 1995) when 
the correction for the disk's contribution is removed. The 
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predicted optical depth at the mean {l,h) = (2.7°, —4.08°) 
of the MACHO fields is ~ 15% lower than in Baade's win- 
dow. Both quoted measurements are consistent with our 
prediction only at the 2 — 2.5c7 level. 

(iii) For 13 clump giants Alcock et al. (1995) find T-e = 
S.QIj ^ (2o"). These values are averages over ~ 12 square 
degrees centered at {l,h) = (2.55°, —3.64°). Our predicted 
optical depth at this position is ~ 7% higher than the result 
for Baade's window. Thus again our estimate is consistent 
with the Alcock et al. result only at the 2<t level. 

A more accurate comparison with the observations 
would involve averaging the optical depth in Fig. 7 over the 
various observed fields, taking into account the number of 
stars monitored in each field, the total time monitored per 
field, etc. 

5 DISCUSSION AND CONCLUSIONS 

We have shown that models of the Galactic bulge/bar and 
inner disk which (i) agree with the COBE/DIRBE L-band 
photometric data, and (ii) have a constant L-band mass-to- 
light ratio Tl, fail to reproduce the reported microlensing 
optical depths to Baade's window by a factor of 2 — 3 the 
predicted optical depths lie 2 — 2.5a below the measured val- 
ues. In this Section, we summarize the steps that lead to this 
conclusion and discuss possible ways out of this discrepancy. 

5.1 The Galactic bar is short, flat, not very 
elongated, and not very close to end-on 

The asymmetry maps of simple analytical bars are domi- 
nated by twin peaks. The location of these peaks varies 
with the orientation, elongation and cixis ratios of the bar in 
a way that is readily understood. In particular, fatter bars 
have wider peaks, and lengthening a bar moves its peaks 
towards higher / and |&|. A near end-on bar can have its 
peaks at similar locations to those of the Galaxy, but then 
its peaks are elongated differently from the Galaxy's. 

Lucy-Richardson iterations that start from one of these 
simple bars considerably improve the fit of the model's 
asymmetry map to that of the Galaxy by deforming the bar 
into an approximately standard shape. As was discussed in 
Paper III, this has axis ratios 1:0.6:0.4. The length of the 
three-dimensional part of this bar/bulge is around 2kpc. 
The part of the bar that extends to ~ 3 — 3.5 kpc is strongly 
flattened to the plane. 

If the Sun-centre line is inclined at (j)o — 25° to the long 
axis of the bar, this standard bar provides a good fit to the 
asymmetry map of the COBE data. But if the bar is oriented 
much more nearly end-on, it provides an unacceptable fit to 
the observed asymmetry map, mainly because the model 
peaks have the wrong shape. 

Moreover, independently of the seriousness with which 
one views the failure of near end-on models to fit the ob- 
served asymmetry map, there can be little doubt as to the 
shape of the near-IR luminosity density in the inner Galajcy, 
because this depends so insensitively on ^o- 

5.2 Bulge mass and microlensing optical depth 
for constant Tl 

If we assume that the mass-to-light ratio Tz, is some fixed 
number, we can determine this number by comparing sim- 
ulated terminal-velocity plots with observed ones. We feel 



some confidence in this calibration because it yields an ex- 
cellent fit to the Galactic terminal velocity curve (Fig. 5) 
and because the simulated {I, v) plots are similar to the ob- 
served {I, v) diagrams in that they show features such as the 
3kpc-arm (Englmaier & Gerhard 1997). 

Given the assumption of constant mass-to-light ratio 
Tl, the mass distribution of the Galactic bar and inner disk 
is rather precisely constrained, and with it the optical depth 
to microlensing t. Nearly independent of the bar angle (/"o, 
T-6 = t/10~^ = 0.8 — 0.9 for main sequence sources and 
r_6 ~ 1.2 — 1.3 for clump giants, where the range of values 
given reflects the uncertainty due to differences between var- 
ious allowed bar models. For a bulge star at 8 kpc distance 
in Baade's window, foreground 'bulge' lenses with galacto- 
centric distances — 3 kpc contribute ~ 55% to the total 
optical depth. These numerical values are for a tangential 
velocity of the LSR of vo = 200 km s"^ if vo = 220kms"\ 
the terminal velocity curve is not fit as well, but for this case 
the inferred mass-to-light ratio and optical depths are 10% 
larger. 

These results arc inconsistent with the measured optical 
depths at the 2a level, and are in conflict with the optical 
depths predicted by some earlier photometric bar models 
- for Baade's window Zhao, Rich & Spergel (1996) predict 
T_6 = 1.1 for the bar only, while Zhao & Mao (1996) predict 
optical depths up to t_6 = 3. 

There appear to be two reasons why our predictions lie 
below those of other photometric models: 

(i) As discussed above, our model luminosity distri- 
bution differs from those used by other authors in three 
important respects: it is more centrally concentrated, 
more strongly flattened towards the plane and less elon- 
gated. These difi^erences reflect the facts (i) that as a non- 
parametric model it fits that data better than do, for exam- 
ple, the widely used Dwek et al (1995) parametric models, 
and (ii) it is based on a much more sophisticated model of 
dust absorption than the foreground-screen model of Arendt 
et al (1994). The central concentration, the large flattening 
and the moderate elongation of our model all conspire to 
diminish, for a given bulge mass, the predicted value of r in 
Baade's window because this field lies 4° from the plane. 

(ii) We determine the mass-to-light ratio Tz, from gas 
velocities rather than from stellar kinematics. In our view 
gas velocities provide the securer normalization because it is 
easier to model the relevant gas dynamics than the required 
stellar dynamics. In particular, it is at present clear neither 
what the line-of-sight distribution of any of the observed 
stellar samples is, nor how anisotropic the velocity ellipsoids 
of these stars are at various points along the line of sight 
(e.g. Sadler, Rich & Terndrup 1996). By contrast, the gas 
streamlines beyond R > 3 kpc deviate from circular orbits 
by only about 5%. 

These differences in the way that we have constructed 
our bulge model lead to our bulge mass being relatively 
small. Integrating over the total NIR luminosity distribu- 
tion, multiplying by the normalization of eqs. (3), and sub- 
tracting the mass of the initial disk models used for the it- 
erations in Paper III yields a bulge mass of Mh — 7.2 — 8.6 x 
lO'' M0 inside 2.4 kpc. By contrast, Zhao, Rich & Spergel 
(1996) derive Mb ~ 2 x lO"' M©, while Zhao & Mao (1996) 
assume 1.8 x 10^° < Mb/M© < 2.8 x 10^°. These masses 
are more in accord with the total bulge plus disk mass in 
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our model, which is Mtot(< 2.4 kpc) = 1.9 x 10^° M©. 

5.3 How to resolve the discrepancy? 

We believe that our models of the inner Galaxy are the most 
carefully constructed models available, and that the lensing 
optical depths that they predict should be treated with re- 
spect. Their ability to give a good qualitative account of 
observed {l,v) diagrams (Englmaier & Gerhard 1997) gives 
one considerable confidence in their fundamental correct- 
ness. Wo trust that this confidence will soon bo tested by 
an examination of their ability to account equally well for 
the data on stellar velocity dispersions in the inner Galaxy 
- complete dynamical models are in an advanced state of 
preparation. But even now it is puzzling that there is such 
significant confiict between the optical depth obtained from 
the microlensing experiments and that obtained from our 
modelling of the COBE and HI terminal velocity data. Sev- 
eral possible causes of this discrepancy stand out. 

The first is that the corrections of Spcrgcl et al (1996) 
for dust absorption are incorrect. As we have emphasized, 
the non-axisymmetric shape of the central luminosity distri- 
bution depends upon asymmetries of < 0.4 magnitudes in 
the adopted dust-free brightness distribution. Consequently, 
small errors in the absorption corrections that have been ap- 
plied to the COBE data could have significant effects on the 
final model, especially in the vicinity of the Galactic plane. 
However, against this possibility one must count the rig- 
orous checks which Spcrgcl ct al showed their corrections 
could pass. The remaining scatter in the dereddened K — L 
colours is only 0.076 mag, which they take as a reasonable 
estimate of the uncertainty in L. Any errors in the ab- 
sorption corrections are likely to be significantly smaller in 
Baadc's window than in the Galactic disk at lower latitudes. 
Notice also that any dumpiness of the dust distribution is 
automatically taken into account in the scaling to the NIR 
reddening. Finally, the experiment reported below shows 
that if the mass in the Galactic plane has been overesti- 
mated, the indirect effects on the optical depth cannot be 
very large. 

Another possible source of error might be the assump- 
tion of eight-fold (triaxial) symmetry assumed in the de- 
projection of the bulge. If the main part of the bulge were 
orientated more end-on than the part that causes the most 
prominent features in the difference maps, the optical depths 
obtained from our models might underestimate the real op- 
tical depths. There is a limit to the amount by which we 
can stretch bulge light along the line-of-sight, however: the 
bulge/bar proper must lie within its corotation radius, which 
from the gas dynamics we know to be ~ 3 kpc. Therefore, 
we can obtain an upper limit on the magnitude of the ef- 
fect of rearranging part of the bulge in the following way. 
We redistribute the surface density in Baade's window that 
arises from the density at line-of-sight distances in the range 
[—3 kpc, 3 kpc] from the center (75% of the total), homoge- 
neously in that distance range. This increases the number 
of sources at the far end and the number of lenses at the 
near end of the bar, so that the total optical depth increases 
by 20% (for Model 2). This is a substantial overestimate of 
the effect of redistributing some of the bulge light, because 
any bulge, whatever its orientation, will be inhomogeneous. 
Thus the error introduced by our assumption of eight-fold 



symmetry must be below 10%. 

The third possibility is that searches for microlensing 
events axe more efficient than their practitioners estimate. 
It is obviously hard for us to comment on this possibility. 
Blending effects have been suggested as a possible cause of 
overestimating the detection efficiencies (Alard 1996) at per- 
haps a ~ 30% level for main sequence stars. However, these 
effects are expected to be less important for clump giants. 

The final possibility is that the L-band mass-to-light 
ratio varies significantly with position in the inner Galaxy. 
Specifically, if T_t were higher above the plane than in the 
plane, the map of optical depth shown in Fig. 7 would have 
too steep a gradient away from the plane. This might be 
the case if the L-band light distribution had a significant 
contribution from young supergiant stars unrelated to the 
bulk of the old stellar population, or if there were significant 
contribution to the L-band light from PAH 3.3 /im or dust 
emission, or if the distribution of stellar mass were more con- 
centrated to the plane than the distribution of gravitating 
mass. 

We have tested this by considering mass models in 
which Tl varies as T_l = Tlo[1 + tanh(|2|/150 pc)]. This 
doubles the mass-to- light ratio at large \z\ with respect to 
the value T^o in the plane. On recomputing the potential 
and choosing T^o so as to keep constant the circular speed 
at 3 kpc, we find that Tlo is 33% lower than the value of 
T L that is given in eq. (3a). Thus the constraints from the 
terminal-velocity curve ensure that the mass of the upper 
bulge increases only by a factor of 1.34 when we assume 
that T is there twice as great as in the plane. The resulting 
values for the optical depth in Baade's window, r_6 = 1-08 
for /3 = — 1 main sequence stars and t-6 = 1.5 for /3 = 
clump giants, still lie well below the measured values. 

Similarly, even models that include 'dark disks' are un- 
likely to give much higher optical depths to Baade's window 
than those calculated here, for the following reason. The 
good fit of our model terminal velocity curve to the observed 
terminal velocities in Fig. 5 argues that the radial distribu- 
tion of mass in the model is approximately correct. Thus the 
radial mass distribution of any additional component has to 
be similar to that of the NIR light, and our derived mass-to- 
light ratio automatically includes its contribution. The main 
remaining freedom is to shift mass out of the Galactic plane 
while keeping the integrated surface density constant. This 
would be most effective about half-way between the Sun and 
the Galactic centre, where the NIR disk has a vertical ex- 
ponential scale-length of hz ~ 150 pc (Paper HI, Fig. 9). 
The line-of-sight to Baade's window passes about ~ 300 pc 
above the plane at galactocentric radius 4 kpc. If the ver- 
tical distribution of mass is also exponential, but with a 
longer scale-length and correspondingly reduced mass den- 
sity in the plane, then the density at height 300 pc can be 
increased by up to a factor 1.36. This optimal value occurs 
for hz — 300 pc. Since thickening the disk reduces the con- 
tribution from lenses near the Sun and has little effect in the 
bulge, we conclude that by thickening the mass distribution 
one cannot increase the optical depth to Baade's window by 
more than ~ 20%. 
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5.4 Conclusion 

Constant mass-to-light ratio models of the inner Galaxy, 
which are consistent with the COBE/DIRBE near-infrared 
photometry and the HI and CO terminal velocity curve, re- 
sult in low optical depths for bulge microlensing, almost in- 
dependent of the orientation of the bar: r_6 = t ~ 0.9 
for main sequence sources, and t-q ~ 1.3 for clump giants. 
These values are inconsistent at the 2 — 2.5cr level with the 
optical depths inferred from microlensing experiments. We 
have discussed several possible uncertainties, including pos- 
sible variations of the near-infrared mass-to-light ratio, but 
none of them appears to be large enough to explain the dis- 
crepancy. 
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